Spin-orbit coupling and perpendicular Zeeman field for fermionic cold atoms: 
Observation of the intrinsic anomalous Hall effect 
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We propose a scheme for generating Rashba spin-orbit coupling and perpendicular Zeeman field 
simultaneously for cold fermionic atoms in a harmonic trap through the coupling between atoms and 
laser fields. The realization of Rashba spin-orbit coupling and perpendicular Zeeman field provides 
opportunities for exploring many topological phenomena using cold fermionic atoms. We focus on 
the intrinsic anomalous Hall effect and show that it may be observed through the response of atomic 
density to a rotation of the harmonic trap. 
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Two important ingredients for manipulating electron 
spin dynamics and designing spin devices in spintron- 
ics [l| are spin-orbit coupling and Zeeman field. For in- 
stance, Rashba spin-orbit coupling (RSOC) (assume in 
the xy plane), together with a perpendicular Zeeman field 
(PZF) (along e^), yield a transverse (along e^) topologi- 
cal Hall current with an applied electric field (along e^). 
Such an intrinsic current was proposed to be one feasible 
explanation of the experimentally observed anomalous 
Hall effect (AHE) and spin Hall effect (SHE) in ferro- 
magnetic semiconductors However, the scattering 
of electrons from impurities and defects in the solid, lead- 
ing to extrinsic AHE and SHE, makes the experimental 
observation of intrinsic AHE and SHE very difficult Q. 

Ultra-cold atomic gases experience an environment es- 
sentially free from impurities and defects, and therefore 
provide an ideal platform to emulate many condensed 
matter models or even observe new phenomena. One im- 
portant recent effort along this line is in study of atom- 
tronics that, in analogy to electronics and spintronics, 
aims to realize devices and circuits using cold atoms 5]. 
One natural and important question in atomtronics is 
how to generate effective spin-orbit coupling and Zee- 
man fields. Great progress has been made recently on 
the generation of RSOC by considering the coupling be- 
tween cold atoms and laser fields [fl-H, which leads to a 
series of important applications |10l4l3| . 

However, the direction of t he g enerated Zeeman field in 
these previous schemes 0, [9l4l3| is in the spin-orbit cou- 
pling plane. Such in-plane Zeeman fields cannot open a 
band gap between different energy branches in the en- 
ergy spectrum. The band gap, together with RSOC, is 
the physical origin of many topological phenomena. For 
instance, the band gap is necessary for the observation 
of the intrinsic AHE [2J. It is also the key ingredient of 
the recently broadly discussed schemes on the creation 
of a chiral Px + ipy-wstve superfluid /superconductor from 
an s-wave superfluid/superconductor |T^. [l5| for the ob- 
servation of non-Abelian statistics and topological quan- 
tum computation [16]. In contrast, a band gap can be 
opened in the presence of RSOC and PZF, but a scheme 
for generating them simultaneously for cold atoms is still 
absent. 



In this paper, we propose a scheme to create RSOC 
and PZF simultaneously for a fermionic atomic gas in a 
harmonic trap. The realization of RSOC and PZF may 
open opportunities for the observation of many topolog- 
ical phenomena in cold atoms because of the non-zero 
Berry phase induced by RSOC and PZF. Here we focus 
on one of them: observation of the intrinsic AHE. In solid 
state systems, the AHE has been observed in transport 
experiments for electrons (i.e., measuring charge currents 
or voltages) . Such transport experiments are not suitable 
for cold atoms in a harmonic trap. We find that the time- 
of-flight of cold atoms in the presence of RSOC and PZF 
can only yield a small asymmetry of the atomic density, 
therefore it may not be suitable for observing the intrinsic 
AHE. Instead, we consider the response of atom density 
to an external rotation of the trap, which corresponds 
to an effective magnetic field for atoms. Because the 
intrinsic AHE is not an quantized effect, the Streda for- 
mula [l7| that was proposed for studying quantum Hall 
effects in cold atoms ^8|] does not apply. We find that 
the atomic density response to the rotation contains not 
only contributions from the anomalous Hall conductivity, 
but also a new term from the orbit magnetic moments 
of atoms that is absent in previous literature [19]. Both 
contributions originate from the topological properties of 
RSOC and PZF. 

Consider ultra-cold Fermi atoms with a tripod elec- 
tronic level scheme (Fig. [T^). States Jl), 12), ]3) are 
three hyperfine ground states, and state ]4) is an excited 
state. The fermi gas is confined in a quasi-two dimen- 
sional (xy-plane) harmonic trap. Along the z direction, 
the atomic dynamics is "frozen" by a deep optical lattice, 
leading to a multiple layered system. The ground states 
]1), ]2), ]3) are coupled with state ]4) by three lasers with 
corresponding Rabi frequencies ftai, ^a2, and ilas- In 
the interaction representation, the single-particle Hamil- 
tonian is 



H = p'^/2m + Vext +Hi- n, 



(1) 



where Hi = /lA 14) (4] - ?i(f7,i 14) (1] + ^i^a K) (2] + 
14) (3] -I- H.c.) describes the laser-atom interaction, 
where A is the detuning to state 14). Vext is the exter- 
nal potential that includes the harmonic trap as well as 
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Figure 1. (Color online) Schematic representation of the light- 
atom interaction for the generation of the effective Hamilto- 
nian ((l}. a, b, c are three sets of lasers. Solid: fiai; Dotted: 
flbi; Dashed: Qci- A is the laser detuning of laser set a, 5 is & 
small shift of the detuning for laser sets b and c from A. (b) 
The configuration of the laser beams. All lasers are uniform 
plane waves. 



Ac = A - 5, with (5 ~ 27r X 80 MHz< A. The 
small shifts of the detunings do not change the effec- 
tive Rabi coupling between different hyperfine ground 
states, but remove the interference among different sets 
of lasers. The optical potentials generated by the laser 
sets b and c are taken as external potentials. With 



suitably chosen Rabi frequencies: 51 
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\/2flQe^^'-\ the Hamiltonian (jSj reduces to Hz 



^b2 

iho {\D2) {Di\ - \Di) (D2I) with ho = Shfll/A. Here 
is the magnitude of the Rabi frequency of the laser b2. 
Within the dark state basis \Di), Hz = hoCTy. Un- 
der a new dark state basis |xi) — {\Di) — i [£'2)) /V^, 
1X2) — \Di) + I-D2)) / \/2, which corresponds to a uni- 
tary rotation of \Di), the Hamiltonian ^ becomes 

He = 7fc^ - ^ + a {k^cFy - kya^) + /iqctz + V" (r) , (4) 



potentials created by other laser fields, ^ is the chemical 
potential. 

The Rabi frequencies fiai can be parameterized as 
rial — sin cos lyse^'^i , Vla2 = i7a sin sin 9?e''^^ , 
17,3 = ^^acos0e'^^ and 17, = 

The diagonalization of Hi yields two degenerate dark 
states: \Di) = sin^se^'^^i |1) - cost/je*'^'^^ |2), jZ^a) = 
cos 6* cos i^e'-^^i |1) -I- cos 6* sin (pe**'32 |2) - sin6'|3), with 
Sij = Si — Sj. We choose a laser configuration il- 
lustrated in Fig. [TJd with the parameters ip = 7r/4, 
9 ~ arctanA/2, 5*.^ = • R, qi = qa (cos ^e^ + sin f e^^), 
q2 = go (cosfe^r - sin|ey), qa = qie^, R = {x,y,z). 
The lasers fiai, ^a2 are in the xy plane, ilaS is along the 
z direction, and \^lai\ = 1 17,21 = |17q3|. Note that these 
three lasers are uniform plane waves, which are different 
from the optical lattices used in a previous scheme [l3| 
to generate RSOC. The effective low-energy Hamiltonian 
is obtained by projecting the Hamiltonian ^ onto the 
subspace of the degenerate dark states spanned by \D.i) 



He = 7^:^ - fi- a [kxOz + kyGx) + V, 



(2) 



where 7 = /2m, a = 79o/\/3, Vij = {Di\Vext\Dj) 
is the effective external potential. The external har- 
monic trap Vtrap = V (r) (^1 is chosen to be spin- 
independent to avoid heating of atoms, where V (r) — 
maj^r^/2, r = (x, y) is the coordinate in the xy plane, ujt 
is the trapping frequency. 

Two additional laser sets b, c with Rabi frequencies Vlbi, 
riei in Fig. [1] induce two Raman transitions between dif- 
ferent hyperfine ground states, yielding an effective cou- 
pling interaction 
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for atoms at the hyperfine ground states 
where F 1 = IULI + l^c2L 2 
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17*2 17c3- We choose the detunings A;,, A^, 



where the third term is the RSOC, the fourth term is the 
PZF. All eight lasers used for the generation of RSOC 
and PZF are uniform plane waves, therefore they do not 
lead to spatial periodic modulation of the atomic den- 
sity. In addition, these lasers propagate only along three 
different directions (the same as other tripod schemes 
d, [13, [m E^D) therefore our scheme should be feasible 
in experiments. 

Under the local density approximation with the local 
chemical potential n (r) = ^ — the Hamiltonian 

@ has two eigenenergies £k± — 7^^ ~ M (r) ± £^0 with 
Eo = a/Ziq + a^A;^. There is an energy gap Eg = 2ho 
opening between two spin orbit bands at fc = (Fig. [2^). 
The intrinsic AHE is nonzero only when the chemical 
potential /i (r) lies inside the gap. 

The dynamics of cold fermi atoms are described by the 
semiclassical equations of motion [20j 



f = aek/5k - k X r^, k = F/n, 



(5) 



for the two sets (b and c) of lasers as A;, 



where F (k) — a'^ho/2 [a'^k^ + HqY ^ is the Berry cur- 
vature in the momentum space and the physical origin of 
many topological phenomena. When there is a nonzero 
external force F along the x direction, an anomalous ve- 
locity -vah — — k X r of atoms is induced along the y 
direction, vah is the physical origin of the intrinsic AHE 
in electronic systems. However, it is difficult to perform a 
similar transport measurement commonly used for elec- 
tronic systems for cold atoms in a harmonic trap. In the 
schemes for the observation of SHE in cold atoms 0- 
the time-of- flight image has been proposed. Because 
of the anomalous velocity vah^ we expect the time-of- 
flight expansion of cold atoms with RSOC, PZF and the 
gravitational force F — mg (along the x direction) has a 
transverse shift along the y direction. 

We assume that the harmonic trap in the xy plane 
1), |2), |3), is suddenly turned off at t = 0, and atoms start to ex- 
Fa = pand, following the equations of motion ([5]). The column 
density along y direction after the expansion depends on 
(5, the initial atom distribution in both momentum and real 
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Figure 2. (Color online) (a) The spin-orbit band structure. 

use the wavelength 

A = 767 nm of the "^^K D2 transition line, the wavevec- 
tor fcfl = 27r/A, the recoil frequency lor — hk%/2m — 
2iT X 8.5KHz, the recoil energy Er = hwR, as the units of 
length, wavevector, frequency, and energy, respectively. Er, 
corresponds to a temperature T = 2fiK. (b) The time of flight 
image of atoms in the presence of PSOC, PZF, and gravita- 
tional force. The unit of the density is kR = 6.7 x 10^ cm~^. 
Ai = 0, T = 20nK, ut = 2tt x IQHz. Solid: £ = 0; Dotted: 
t = 19 ms. 



spaces 



\Vf.T,t) = J d?v, J 



(2^ 



Siyf-y,^Ayit))f, (6) 



where T is the temperature, t is the time of flight, 
(jji) is the initial position of atoms, and /(k, ri,T) ~ 
1/ (exp [(eic — /i (r^)) /ksT] + 1) is the Fermi-Dirac dis- 
tribution of atoms, yf is the position of atoms at time t. 

^yit) = /':^^'^'<(Tfe+r.(fc;,fc,)) is the dis- 
tance of flight. 6 (y) is the delta function. 

In Fig. [2}d, we plot the column density n{yf,T,t) at 
two different times. At t = 0, the initial column density 
is symmetric along the y axis. After a time t, n{yf,T, t) 
becomes asymmetric because of the anomalous velocity 
of atoms. However, the expansion dynamics of atoms 
is dominated by the first term in Ay (t), and the differ- 
ences of the atomic densities at ±yf correspond to only a 
small percentage (below 3%) of the total density. There- 
fore it may be hard to observe them in experiments. The 
difference between electrons and cold atoms comes from 
the fact that de]s_/dky in an electronic system does not 
contribute to the overall transverse motion (the observed 
Hall current is purely from the anomalous velocity), while 
the initial atom velocities d£\^/dky dominate the expan- 
sion process and the asymmetry of the column density 
should be small in an atomic system. Therefore we need 
develop other techniques for the observation of the in- 
trinsic AHE. 

Recently, the response of atomic density to a rotation 
of the trap has been proposed for measuring quantum 
Hall conductivit y fo r cold atoms based on the well-known 
Streda formula cr — dn/dB. Here the rotation for 
atoms is equivalent to the magnetic field for electrons. 
However, the Streda formula does not apply to the AHE 
because it is not quantized. Nevertheless, the density 
response to the rotation still contains rich information 
about the intrinsic AHE. Consider a rotation of the har- 
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Figure 3. (Color online) Plot of the atom density n (ro, /i, T) 
with respect to km and /x. The units are the same as that 
in Fig. El /lo = 0.05-Efl. = I/ZEr. Sohd and dotted: 

T = 2nK; Dashed and dash dotted: T = 200717^. (a) Solid 
and dashed: ^ = Q; Dotted and dash dotted: ^ = — 0.01-Eii. 
(b) Solid and dashed: = 0; Dotted and dash dotted: vj = 
O.OSwfl. 



monic trap along the z axis [21| , the Hamiltonian can be 
written as 

H = h^q^ /2m + a [q^iyy - qy(Jx) 

+/iocr^ + TO (cj^ _ tn^) rV2 - (7) 

in the rotation frame, where = Kk — mwz x r is the 
mechanical momentum of atoms [ioj, w is the rotation 
frequency of the trap. The density of atoms is 

n (r) = {2Try^ J d\ [1 -I- mwT./H] f (q, r,ro) , (8) 

where mzuTz/h is a correction to the well-known con- 
stant density of states 1 / {2tt) in the presence of 
nonzero Berry curvature fields and the rotation \x§l,. 
f (q, r,m) = 1/ (exp [{Sm ~ M (r)) /ksT] + 1) is the Fermi 
Dirac distribution of atoms, where the energy of atoms 
£m = £q ^ MzVJ contains a correction —M^vj — 



hi) 



known as the magnetization 



energy for electrons in the solid state [2C|] . 

We assume the rotation frequency w of the trap is 
slowly increased to keep the same temperature of the 
system. The local chemical potential can be fixed by 
increasing ujt with vj. The response of the atom density 
to VJ is 



dn 



{2nf 



mnjT,\ df 



(9) 



for fixed fi and T. At T = and = 0, Eq. 
to 



dn 
dm 



TO 



(2^ 



1 



47r dsq/dq'^ 



reduces 



(10) 



Here the first term is the anomalous Hall conductivity 
axy for cold atoms. In the parameter region |/i(r)| < ho 
(i.e., the chemical potential lies in the band gap), it 

yields a^y = 



1 ^ hi^l ^oP-q^p (r) -f /iq /47r?i, where 
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Figure 4. (Color online) Plot of dn (vj, \i, T) /dzu with respect 
to zu and fi. The unit of dn {lu, fj,,T) /dzu is taken to be 
m/2Tih. The other units are the same as that in Fig. (2] 
h(i = 0.05Sfl and a = (a) T = 2 nK (b) T = 200 

nK. 



the Fermi wavevector qp {r) is obtained from ec^p — ^J■ir)■ 
In the parameter region 01^/7 » ho, qp ~ Il^ and 
(jxy « mjA-'Kh. The second term 0-3:^, originating from 
the non-zero orbit magnetic moment M^, is an addi- 
tional contribution to dnjdvj that was missing in the 
previous literature 19] for electron systems. Eq. ([10]) 
is a generalization of Streda formula for the anomalous 
Hall effects. By varying parameters and measuring the 
density response, we can extract information not only 
about the anomalous Hall conductivity, but also the 
magnetic moment that is generally hard to measure in 
solid state systems. In the parameter region ^ /iq, 
a^y ~ ^mhol%'nho? ^ mjA.'Kh « cjxy Therefore a^y 
dominates in Eq. ([9]) in this region, and the density 
response dnjdvj yields a rough measurement for the 
anomalous Hall conductivity. 

We numerically calculate the density n and density re- 
sponse dnjdvj as functions of the parameters (p^za^T), 
and plot them in Figs. [3]and|lfor T = 200 nK and 2 nK. 
The presence of the harmonic trap changes the chemical 
potential at r by — m (w^ — nj^) r^/2. In a realistic ex- 
periment, the effective trapping frequency \J uJl ~ vj^ in 



the presence of rotation may be slightly different from 
the initial trapping frequency Wjo without rotation to 
keep the same temperature of the system (2^ . This can 
be overcome by comparing the densities at different spa- 
tial points r, r* such that /i — \m (ujf — ro^) r*^ = M — 
imwjo^^ ^'^ keep the same local chemical potential [2^ . 
With this method, we can measure the density response 
to the rotation with the fixed temperature and chem- 
ical potential. In addition, the rotation of the system 
requires an asymmetric harmonic trap [2ll |. which does 
not affect our results because it only change the spatial 
positions for the measurement of the density change at a 
fixed chemical potential through a different local chemi- 
cal potential dependence /i (r) = /i — TO[(a;j^ — tu^) -I- 
(wjj^ — zu'^) 2/^]/2. We adopt a set of parameters: a'^/j = 
Er/3 = 27r?ix2.8KHz, Hq = 27rfix425 Hz, ujto = 27rx50 

Hz, tUniax = O-ObLdR — 2tT X 4:25Hz, OJt ~ y^Ld^Q+W^, 

LOt — vj = ujIqI {ujt + w) ^ 2Tr x3Hz. At r = O.lXuJu/TTujto 
— 4.2/im, the chemical potential n changes by O.OlEji. 
From Figs. [3] and HI we see a maximum density change at 
the order of 3 x 10^ cm~^ can be observed with a rotation 
frequency romax, which corresponds to about 10% of the 
total density and can be observed in a realistic experi- 
ment. The density variation at a medium temperature 
T — 200 nK is at the same order as that at a low temper- 
ature T — 2 nK. Note that the multiple layer structure 
induced by the optical lattice confinement along the z 
direction can further enhance the signal. 

In summary, we propose a scheme to create RSOC and 
PZF simultaneously for cold atomic gases. We show that, 
by measuring the atomic density response to a rotation 
of the trap, the intrinsic AHE can be observed for cold 
fermionic atoms in a harmonic trap. We emphasize that 
the creation of RSOC and PZF brings new opportunities 
for studying many topological phenomena, such as chiral 
p-wave superfluids, anomalous and spin Hall insulators, 
etc. 
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